Resistivity of Inhomogeneous Superconducting Wires 
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We study the contribution of quantum phase fluctuations in the superconducting order parameter 
to the low-temperature resistivity p{T) of a dirty and inhomogeneous superconducting wire. In 
particular, we account for random spatial fluctuations of arbitrary size in the wire thickness. For a 
typical wire thickness above the critical value for superconductor-insulator transition, phase-slips 
processes can be treated perturbatively. We use a memory formalism approach, which underlines the 
role played by weak violation of conservation laws in the mechanism for generating finite resistivity. 
Our calculations yield an expression for p{T) which exhibits a smooth crossover from a homogeneous 
to a "granular" limit upon increase of T, controlled by a "granularity parameter" D characterizing 
the size of thickness fluctuations. For extremely small D, we recover the power-law dependence 
p(T) ~ T°' obtained by unbinding of quantum phase-slips. However in the strongly inhomogeneous 
limit, the exponent a is modified and the prefactor is exponentially enhanced. We examine the 
dependence of the exponent a on an external magnetic field applied parallel to the wire. Finally, 
we show that the power-law dependence at low T is consistent with a series of experimental data 
obtained in a variety of long and narrow samples. The values of a extracted from the data, and the 
corresponding field dependence, are consistent with known parameters of the corresponding samples. 



I. INTRODUCTION 

Transport in superconducting systems of reduced di- 
mensions (thin films and wires) is known to be strongly 
affected by fluctuations in the order parameter. A promi- 
nent manifestation of the role of fluctuations is the fi- 
nite electrical resistance of narrow superconducting (SC) 
wires at any finite temperature T below the bulk critical 
temperature Tc, established when the wire thickness d is 
reduced below the superconducting coherence length^ ^. 
The finite voltage drop along the wires is generated by 
phase-slips: these are processes whereby the supercon- 
ducting phase slips by 27r at points on the wire where 
superconductivity is temporarily destroyed. The result- 
ing voltage is related to the rate of phase slips via the 
Josephson relation. The pioneering theoretical studies 
of this phenomenon^ accounted for thermally activated 
phase-slips across barriers separating metastable phase 
configuration states, corresponding to local minima of 
the Landau-Ginzburg free energy. This theory turns out 
to be supported by experimental studies^, in particular 
providing a good fit of the resistivity p{T) slightly be- 
low Tc- As T is lowered further, p{T) is exponentially 
suppressed and becomes practically undetectable. 

The advance of nanostructure fabrication techniques 
during the 1980's opened up the possibility to study nar- 
row wires of smaller diameter, down to a few 100 A. The 
consequent weakening of superconductivity leads to en- 
hancement of the rate of phase-slips, and results in a 
measurable finite resistance even at T <C Tg. In this 
low T regime, however, thermal activation is considerably 
suppressed, and the dominant mechanism for phase-slips 
becomes quantum tunneling. The first experimental in- 
dication of quantum phase-slips (QPS) has been seen in 
thin strips of Indium by Giordanc'^. The curves of p{T) 



vs. T exhibit a 'kink' at some T* smaller than Tc, be- 
low which the decrease of p{T) upon lowering T becomes 
more moderate. For T < T*, the fit to the theory of 
Ref. fails. Motivated by the assumption that quan- 
tum tunneling dominates over thermal activation in this 
regime-, Giordano fitted the data to a phenomenological 
expression in which the temperature fc^T in the thermal 
activation rate ^ ^-AF/ksT replaced by a "character- 
istic frequency" given by h/T^Q, where t^q is the re- 
laxation time in the time-dependent Landau-Ginzburg 
theory. While this phenomenological expression fitted 
the data quite successfully, it is not justified by a rig- 
orous theoretical derivation. In particular, the Landau- 
Ginzburg theory (which is employed in the calculation of 
the barrier hight AF) is based on an expansion in the 
close vicinity of Tc, and is not expected to hold far below 

Tr. 



Subsequent theoretical studies of the QPS contribu- 
tion to resistivity below Tt^i^i^i^ d^did^ have addressed 
the quantum dynamics of phase fluctuations in the low- 
temperature limit (T ^ Tc), yielding a radically differ- 
ent behavior of p{T) vs. T. At low T, the magnitude of 
the SC order parameter field is approximately constant 
and its dynamics is dominated by phase-fluctuations. 
The corresponding T — > effective model is a (1-1-1)- 
dimensional XY-model, in which the topological excita- 
tions are QPS and anti-QPS (vortices and anti-vortices 
in space-time )S.. Upon tuning a stiffness parameter 
(which in particular is proportional to the wire diam- 
eter d) below a critical value, the system undergoes a 
quantum Kosterlitz-Thouless (KT)^'^ transition from a 
SC phase to a metallic phase driven by the unbinding 
of QPS-anti QPS pairs. The latter are responsible for 
a finite resistivity at any finite T even in the SC phase, 
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where 



p(T) - po 



(1) 



with a > 0. This power-law dependence was predicted 
both for a granular wire composed of weakly coupled 
SC grains separated by a tunnel barrier (a Josephson 
junction array)^'^^ and for the case of a homogeneous 
wire in the dirty limit®'^'^^; however, the exponent a and 
the prefactor po are different (the latter, in particular, 
strongly depends on the microscopic details). 

The quantitative estimate of po and the consequent size 
of the resistivity in realistic systems have been discussed 
in detail in the theoretical literature, and led to a de- 
bate among the authors^^ regarding the relevance of the 
homogeneous wire limit to the experimental system of 
Ref. 4]. However, we arc not aware of a direct attempt 
to test the validity of the power-law scaling of p{T) as 
an alternative to Giordano's phenomenological formula. 
More recently, several experimental group s^^'^^i^^'^^i^^ 
have reported QPS-induced resistivity in a variety of SC 
nanowire samples with d of order 100 A and below. The 
data in certain samples also exhibit a transition to a nor- 
mal (weakly insulating) state below a critical wire di- 
ameter dc, consistent with the theoretical expectations. 
In the SC samples, the data was fitted by an effective 
circuit which accounts for both quantum and thermally 
activated phase-slips contributions to p(T). Here as 
well, Giordano's phenomenological expression was imple- 
mented to describe the QPS term in the fitting formula, 
rather than a power-law ~ T" (see, however, the un- 
published notes in Ref. [l^ and 0). The comparison 
between theory and experiment is therefore not entirely 
settled. 

In the present paper, we introduce a derivation of 
QPS-induced resistivity which enables a systematic ac- 
count of all possible scenarios in a realistic experimental 
setup, and carry out a comparison with the experimen- 
tal data. We focus our attention on dirty SC wires, and 
account for non-uniformity of arbitrary size in the wire 
parameters, in particular the diameter d. Our approach 
implements the memory function formalism, which di- 
rectly relates the mechanism responsible for generating 
finite resistivity to the violation of conservation laws of 
the unperturbed Hamiltonian Hq (see below) , describing 
non-singular fluctuations. The formalism also underlines 
the interplay of phase-slip processes and disorder. We 
consider two distinct types of disorder: the first corre- 
sponds to impurities in the underlying normal electronic 
state, characterized by a mean free path £ <^ ^, and the 
second type is associated with the random spatial fluctu- 
ations in the wire diameter on a length scale of order ^, 
whose size is characterized by a "granularity parameter" 
D. We find that p{T) exhibits a smooth crossover from a 
homogeneous to a "granular" behavior upon increase of 
T. The latter is expected to dominate in most of the mea- 
surable range of T due to an exponential enhancement of 
the prefactor po in Eq. ([T]). In addition, we examine 



the dependence of the exponent a on an external mag- 
netic field applied parallel to the wire. The model and 
the main steps of our calculation are described in Sec. 
im details of our derivation of the resistivity within the 
memory approach are given in Appendix A. In Sec. IIIII 
we show that the power-law Eq. ([1]) is consistent with 
the experimental data obtained in a variety of samples, 
provided the wire is sufficiently long. The values of a 
extracted from the data, and the corresponding field de- 
pendence, are consistent with known parameters of the 
corresponding samples. Our conclusions are summarized 
in Sec. [IVl 



II. THE MODEL AND DERIVATION OF 
PRINCIPAL RESULTS 

We consider a long and narrow superconducting (SC) 
wire of length L and cross section s = TT{d/2)'^, such that 
L ^ and d <C £,,Xl where ^ is the superconducting 
coherence length and Xl the London penetration depth. 
Fluctuations in the SC order parameter are therefore ef- 
fectively one-dimensional. As a first stage we assume 
that the wire is homogeneous, and the SC material is in 
the dirty limit where the mean free path in the under- 
lying electronic system obeys i An effective model 
in terms of the SC order parameter field is obtained as 
a result of integrating over the electron fields^ii^. At low 
temperatures T <^ Tc (where Tc is the bulk SC tran- 
sition temperature), we assume that fluctuations in the 
magnitude of the order parameter are suppressed, while 
the quantum dynamics of phase fluctuations (j){x) is de- 
scribed by the Sine-Gordon Hamiltonian® 



where 



Hn = - 



jjps ^ 



L/2 
-L/2 



dx 



c 



H^ 



L/2 



L/2 



dx cos(2n6 



(2) 



(3) 



(4) 



Here and throughout the rest of the paper we use units 
where ?i = fc^ = 1; H(x) is the field conjugate to (t){x), 
satisfying [4>(x), H(a;')] = i5{x' — x), which physically rep- 
resent fluctuations in the number of Cooper pairs and is 
related to the fleld via dxO = ttH. C is the effective ca- 
pacitance per unit length, Ug is the (three-dimensional) 
superfluid density, and e, m are the electron charge and 
mass, respectively. In Hf^^ , y = ex.p{— Score} is the fu- 
gacity of quantum phase-slip (QPS), where Score is the 
action associated with the creation of a single QPS (of 
winding number 1) due to the suppression of the SC or- 
der parameter in its core; the characteristic velocity v is 
given by 
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n is the winding number ("charge") of a QPS, and a is a 
short distance cutoff ~ ^. Note that the Hamihonian Eq. 
([2|) is the same effective model as in Ref. Q expressed in 
a dual representation^. 

The Hamiltonian Hq describes the non-singular phase 
fluctuations, which are characterized by a free mode (the 
Mooij-Schon mode^^) propagating at a velocity v. It can 
be recast in the familiar Luttinger form 



Ho = V 



dx 
2^ 



where v is defined in Eq. ([5]) and 



K = 



TT \snsC 



1 



1/2 



(9.0)' 



(6) 



(7) 



Note that we also assume the wire to be sufficiently long 
so that L is large compared to f/T, in which case it 
can be practically taken to be infinite. The nature of 
the T — Q fixed point of H [Eq. ^] depends crucially 
on the Luttinger parameter K. In particular, as noted 
by Ref. ^ , the system undergoes a quantum Kosterlitz- 
Thouless (KT)^^ transition from a SC phase at K < 
Kc {Kc ~ 1/2) to a metallic phase &t K > Kc- For 
given material parameters, the transition can be tuned 
by varying the wire cross section s below a critical value 
Sc, related to through Eq. ([7]). We hereon focus on the 
SC phase corresponding to s > Sc, in which all the terms 
HPf [Eq. (dl] are irrelevant and the T = fixed point 
Hamiltonian is Hq. Indeed, it describes a true SC state, 
where the resistivity p(T) vanishes in the limit T — > 0. 

A key feature of the SC state is that the charge current, 
dominated by the superconducting component 



Je. 



m 



^ a A 2e 
dxd^cj) = — — 

TT K 



dxdx4i , 



(8) 



is an almost conserved quantity. Formally, this is mani- 
fested by the vanishing of its commutator with the low- 
energy Hamiltonian Hq, which implies that in the ab- 
sence of the phase-slip terms H'p^ , the current cannot 
degrade {dtJe = 0) and hence the resistivity vanishes. 
The leading contribution to p{T) can therefore be ob- 
tained perturbatively in H^'^ . As we show in detail in 
Appendix A, the calculation of p{T) can be viewed as a 
particularly simple application of the more general mem- 
ory matrix approac h^'^'^^i^^'^'^ , which directly implements 
this insight. In essence, this approach incorporates a re- 
casting of the standard Kubo formula for the conduc- 
tivity matrix (a highly singular entity in the case of an 
almost perfectly conducting system) in terms of an ob- 
ject named a "memory matrix" . The latter corresponds 
to a matrix of decay-rates of the slowest modes in the 
system, and is perturbative in the irrelevant terms in the 
Hamiltonian, in particular all processes responsible for 
degrading the currents and hence generating a finite re- 
sistivity. The separation, in this approach, of the slow 
modes generated by the irrelevant operators around Hq 



from the fast modes allows a controlled approximation as 
the temperature is lowered and provides a lower bound 
on the conductivity^!. 

Our derivation of the d.c. electric resistivity (see Ap- 
pendix A for details) yields the following expression: 



P{T) = , 

ee 

where Xee is the static susceptibility 

Xee = 7^{Je\Je) 



(9) 



(10) 



and Mee (the memory function) can be expressed, to 
leading order in the perturbations H'p^ , in terms of cor- 
relators of the 'force' operators 

F;,,„=z[Je,i/r] , (11) 

which dictate the relaxation rate of the current Je via 

dtJe^Y^^kn- (12) 

n 

This yields 

Af,,(r) « Y.Kln. (13) 

n 

( J?^ • fP^ \0 ( • \^ 

\ ps^ni ps,n/oj \ ps^ni ps.n/ uj—O 



lUJ 



in which {F"^ ] F'^Yl^ is the retarded correlation function, 
the expectation value being evaluated with respect to Hq. 

To leading order in perturbation theory, we evaluate 
the expectation value in Eq. (fTO|) as well with respect to 
the low energy Hamiltonian iJp- This yields 



1^ 



(14) 



Using Eqs. ([8]), ([4]) and (fTTj) we obtain an expression for 
the force operator 



Ks,n = ^ ;^ 2^" / ^^^^"^ (2"^ 



(15) 



Inserting in Eq. (|T3|) . we find 

4i(4nei;2y»)2 



X I dx dttlm{Cpt ,,{x,t)} , (16) 



where 



Ge{x,t) = (0(a;,i)0(O,O))° 



(17) 
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in which the Green's function Gg{x^ t) at finite T is given 



G8{x,t) = 



AK 



iraT /v 



sinh{7rr(a; — vt -\- ia)/v} 
naT/v 



(18) 



ih{7rr(a 



vt 



i)/v} 



To find the leading T-dependence of p{T) for small y and 
{aT/v), we neglect the contributions of n > 1. Substitut- 
ing the resulting M^e combined with Xee from Eq. ^ 
in Eq. ([9]), we obtain 



^ ^ (2e)2a T^l/K) ^ 



2naT 



-1 2/i<'-3 



V 



(19) 



where r(a;) is the Gamma function. This essentially re- 
covers (up to a numerical prefactor) the result of Ref. Q , 
which indeed corresponds to the homogeneous wire limit. 
Note that the resulting p{T) exhibits a SC behavior as 
long as K < 2/3, in accord with the renormalization 
group analysis of Ref. [13] ■ 

We next turn our attention to the more realistic situ- 
ation, where inhomogeneities along the SC wire are al- 
lowed. Random fluctuations are possible in all the wire 
parameters, and in particular the diameter may varry in 
space leading to a local cross section s{x), which can be 
assumed to be a random function of x. The most promi- 
nent modification of the Hamiltonian in the presence of 
such spatial fluctuations is manifested in the fugacity y, 
which depends exponentially on the wire diameter via 
the core action Score- Consequently, it becomes space- 
dependent, i.e. y = y{x) where 

y{x) = exp {~Score{x)) = cxp (-5*0 - A(x)) (20) 

where A(a;) is a random correction to the uniform core 
action Sq. The leading (n = 1) phase-slip contribution 
to the Hamiltonian [Eq. ^] now becomes 



Hps — 



-2v 



dx y{x) cos(20) 

2 . exp(-A(a;))cos(26') , (21) 

where yo = exp(— So). Assuming a Gaussian distribution 
of the random function A(a;) 



a 

-2yov 
a 



P[A] = AAcxp<^ - 



2aD 



dxK^{x) 



(22) 



yields the disorder averages A(a;) = 0, A(a;)A*(a::') = 
aD5(x — x') and 



A(x)g-A(x') = / I?AF[A]e^(")e^("') 



+ a6{x — x') [e 



j\ r„2D 



(23) 



The last result is derived from the discrete version of 
the above functional integral, which leads to Sue^^ + 



(1 — 5ij) [here 5ij = a5{xi — Xj)]. Note that the pa- 
rameter D characterizes the degree of granularity in the 
SC wire, with ^/D proportional to the typical amplitude 
of spatial fluctuations in the wire cross section. 

The inhomogeneous phase-slip term Eq. (|2T|) modifies 
the expression for the force operator 



2 /> 

Pps^^ri^^ J dxyix) sin (20) 



(24) 



Substituting in Eq. p3p and performing the disorder 
averaging using the correlation function (p3)) we find 



A/' 



MS 



where 



(25) 



M„„ = 7—; e 



ML 



ALjAev^yo)^ 



dx / dttlTsi{Cll^^{x,t)} 

-OO JO 



,2D 



dttlm{C;i,{0,t)} 



and Cp^ i{x, t) is given by Eq. pT]) for n—\. This yields 
the leading T-dependence of the resistivity for arbitrary 
granularity D in the SC wire: 



Ah = 



(2e)2a 



2TiaT 



2/A--3 



(26) 



Ahe 



D 



A„ 



2TiaT 



TT^T\II2K) 



A,^ 



7r2r2(l/X) 
AT{2/K) 



For extremely weak granularity where D <C 1, the first 
term in (I26p dominates, and the homogeneous result Eq. 
([in]) is recovered. However, for D of order 1 or more, the 
second term is exponentially enhanced by the factor e^^, 
yielding 



{2eYa^'y°' 



2TraT 



V 



2/K-2 



(27) 



This approximation is consistent with earlier predictions 
for granular SC wires^iii. Indeed, it indicates that the 
phase-slips dominating the resistivity occur at narrow 
constrictions in the inhomogeneous wire. Our more gen- 
eral expression Eq. ((26)) implies that for a fixed D, the 
resistivity p{T) vs. T exhibits a crossover from a "homo- 
geneous" power-law behavior [Eq. ^ with a = 2/K — 3] 
to a "granular" limit with a higher exponent a = 2/K—2. 
In most realistic systems we expect this crossover to oc- 
cur at very low temperatures: this is in view of the ex- 
ponential dependence on the granularity parameter D. 
The typical core action 5*0 has been estimated to be of 
order 10 or large r^'^i^^ ; this implies that even small ir- 
regularities 10%) in the wire diameter corresponding 
to Z? > 1 are sufficient to enhance the prefactor of the 
second term in (j26p by at least an order of magnitude. 
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III. COMPARISON WITH EXPERIMENTAL 
DATA 

As shown in the previous section, the resistivity of a 
SC wire with a moderately non-uniform diameter [Eq. 
(p7)) ] far enough below is expected to be well approxi- 
mated by a power-law T-dependence of the form llj . In 
particular, the exponential enhancement by the granular- 
ity parameter D partially compensates for the exponen- 
tial suppression by the fugacity y = e~'^=°''», and conse- 
quently the prefactor [po hi Eq. ([1])] becomes comparable 
to the observed resistivity in a number of experiments. 
In this section we test the relevance to available experi- 
mental data by a direct attempt to fit a power-law. In 
the cases where the fit appears to be reasonably good, 
we can extract the Luttinger parameter K from the ex- 
ponent a = 2/ K — 2 and compare it to an independent 
estimate based on the sample parameters. According to 
Eq. ([7]), K depends in particular on the wire diameter 
d ^ ^/s, the capacitance per unit length C and the su- 
perfluid density n^, which can be expressed in terms of 
the London penetration depth A^: 

K^^^^. (28) 

Here A is a numerical factor which depends on the ge- 
ometry of the wire cross section (e.g., A = y/n in Ref. [J| 
and A = 2 in Ref. [17] ) , and C is typically of order unity^ 
and depends very weakly (i.e. logarithmically) on d; we 
hereon regard it as a fitting parameter. The dependence 
on Al implies that a smooth tuning of K is possible by 
application of a magnetic field parallel to the wire. In- 
deed, this was done by Altomare et al. (Ref. [l^), as will 
be discussed in more detail below. 

We have considered data obtained in several different 
experimental setups: Indium strips studied by Giordano 
(Ref. Q), MoGe wires deposited on carbon nanotubes 
studied by the Harvard group (Ref. [TsUlGj ). more recent 
studies of tin (Sn) nanowires (Tian et al., RefJlTf) and 
long Aluminum wires (Altomare et al., Ref. [3). All of 
these indicate substantial deviation from the thermal ac- 
tivation theory^., attributed to QPS. In addition, we note 
that in most of the SC samples involved in those exper- 
iments the normal resistivity is too low to be considered 
in the strictly granular limit (where distinct grains are 
weakly coupled), although the wire diameter is likely to 
be non-uniform with varying degree of nonuniformity. 
The samples of Ref. (Tslllq exhibit a rather complicated 
behavior, following from a number of reasons. First, most 
of the wires fabricated in this particular methods are rela- 
tively short, hence finite size effects interfere with the T- 
dependence, and the dynamics of phase-slips is crucially 
affected by their backscattering from the boundaries'^^. 
In addition, some of the nanotube substrates are not in- 
sulating, and their (unknown) resistance complicates the 
fitting by additional parameters. Hence, although some 
of these samples can be fitted reasonably well by Eq. ([1]), 
a more detailed quantitative analysis has focused on the 



other experimental papers. 
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FIG. 1: Data points obtained from Fig. 1 of Ref. [3| pre- 
senting the resistance vs. T, re-plotted on a log-log scale; 
the solid lines correspond to fitting functions of the form T" . 
(a) In the low-T (QPS) regime: blue circles (red squares) 
correspond to the samples of diameter 410A (505 A); the cor- 
responding exponents are a — 25.2 ± 0.2 and a — 27.85 ± 0.4, 
respectively, (b) The 410A and 505A diameter samples in 
the high-T regime, and all data points of the 720A diame- 
ter sample (green diamonds); a = 89 ± 13, a = 125 ± 8 and 
a = 260 ± 7, respectively. 

We first consider the data of Ref. [J- Fig. ([T|) presents 
the resistance of two samples as a function of T on a 
log-log scale. Note that the length of the shortest wire 
studied is 80 /zm. In comparison, the effective length 
set by the temperature scale Lt = {v/T) for the rele- 
vant T 4K is oi order 10 fim; this follows from the 
estimate v ^ cd/X^, where c is the velocity of light and 
the measured A^ « 1300A obtained for Indium films of 
comparable thickness'^^. The wires therefore fulfill the 
long wire condition L 3> Lt, and can be considered good 
candidates for testing the scaling of p{T) with T. In the 
low T regime [Fig. (HJa))], the fit to a power-law is very 
good'^-. In addition, the values of a corresponding to 
the two samples of diameters 410A and 505A yield the 
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Luttinger parameters K = 0.0735 and K = 0.067, re- 
spectively. These are consistent with the values obtained 
by inserting the relevant d and Al in Eq. ([28| . provided 
C« 2. 

In principle, one may argue that a fit by a power-law 
with a high exponent (especially within a limited range 
of T) is not easily distinguishable from an exponential 
function. However, to put the above analysis to a further 
test, we tried to fit the data obtained in the high T regime 
(in the close vicinity of Tc) to a power-law as well. In 
this regime, there is no justification for such a functional 
dependence on T, as the thermal activation theory for 
phase-slips^ is expected to work much better. Indeed, 
the result depicted in Fig. ((Hb)) indicates an obvious 
failure of a T" trial function: the best fit to a straight line 
on a log-log plot yields an anomalously large exponent (a 
ranging from 89 to 260) , and a systematic deviation from 
the fitting function. This is in sharp contrast with the 
situation in the low-T regime, and therefore strengthen 
our confidence that in the latter case, the analysis is valid 
and provides a meaningful confirmation of the theoretical 
prediction. 
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Log[T(K)] 

FIG. 2: Data points obtained from Fig. 1 of Ref. present- 
ing the resistance vs. T, re-plotted on a log-log scale: sample 
diameters are 20 nm (blue circles), 40 nm (red squares), and 
60 nm (green diamonds) . The solid lines correspond to fitting 
functions of the form T" , with a = 5.38 ± 0.2, a = 5.16 ± 0.2 
and a = 31.2 ± 1.3, respectively. 

We next consider the data of Tian et al (Ref. [l3|)- The 
results obtained from resistivity measurements in three 
different tin (Sn) nanowires of diameters ranging between 
20 nm and 60 nm are depicted in Fig. ([2]). Similarly to 
Giordano's data, in all three samples a "kink" is observed 
at some temperature T* < Tc, below which the decrease 
of i? as T is lowered becomes more moderate. In the 
two thinner samples, the low-T section of the data is 
quite noisy, and R(T) seems to saturate. This behavior 
is possibly due to serial contribution from the contacts. 
Otherwise, however, the power-law dependence appears 



to be a good fit for T <T*. 

To test the validity of Eq. ([^5]) in this system, one 
requires an independent estimate of A^. A direct mea- 
surement is not available, however, an indirect estimate 
of the ratio Al /d based on a measurement of the effec- 
tive critical magnetic field^ yields XL/dK, 1.5 for o? w 70 
nm, which is close in diameter to the thickest sample of 
Fig. ^. This is consistent with Eq. ([28]) for C « 0.98, 
a quite reasonable estimate for the capacitance. Unfor- 
tunately, reliable data on Al /d in the thinner samples is 
not available. A naive extrapolation of the d-dependence 
of Al to lower values of d based on the Landau-Ginzburg 
theory25 (Al ~ l/Vd) yields K = 0.315 for the 20 nm 
sample, a reasonable approximation to the value of K 
obtained from the fit {K = 0.274), however the expected 
scaling with d fails in the case of the 40 nm sample. In- 
deed, as shown explicitly in Ref. 34], the naive scaling 
AL(rf) works well above 60 nm, but breaks down for the 

= 40 nm wire (no data is given for even thinner sam- 
ples) . 

We finally focus our attention on the most recent ex- 
perimental work of Altomare et al. (Ref. [31). This 
group has studied long and thin Aluminum wires, which 
appear to be ideal candidates for comparison with the 
QPS theory. In addition, the application of a magnetic 
field H enables a continuous tuning of the superfluid den- 
sity in a single sample of fixed d. Following Eq. ([7]), we 
therefore expect the ff-dep endence of K to be related to 
Ug via 1/K{H) ~ ns{H). The functional dependence 
of ng{H) on H can be derived from a simple calcula- 
tion using the Landau-Ginzburg theory at finite mag- 
netic field^. This yields 

Ctt2 \ \r\ 

where He is the bulk critical field, and M a numerical 
factor which depends on the geometry: e.g., JV ~ 2\/6 in 
a thin film of thickness d, and Af — 8 in a cylindrical wire 
of diameter d. marks an estimated critical field where 
the inverse Luttinger parameter is expected to vanish 
according to 

1 1 f H'\"' 

The unpublished notes included in Ref. [l^ enable a 
direct comparison of the experimental data to the above 
prediction. The experimental results include the Ohmic 
resistance as well as current -voltage characteristics which 
include a non-linear part. In the notes, it is shown that 
the latter can be fitted to a power-law V ^ I'^ . This is 
suggested by the authors as a plausible alternative to the 
phenomenological exponential expression used as a fitting 
function in the published Letter. Such power-law in the 
I — V characteristic is consistent with a T-dependence 
of the Ohmic resistance of the form R{T) ^ T", with 
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FIG. 3: Data obtained from the unpublished notes (Ref. [l^), 
plotted in terms of the inverse Luttinger parameter {1/K) as 
a function of magnetic field H . All data points correspond to 
the lowest temperature measurements {T = 0.4 K). The solid 
curve corresponds to the formula Eq. (|30p . with _ff|| ~ 1.1 T. 



a = V — 1. The exponent v is then extracted from the 
data for different values of the magnetic field, and plotted 
as a function of H . Using the linear relation of v to 
I/if, we plot the data in the form depicted in Fig. ([3]), 
where it is fitted to a trial function of the form pO|) . The 
fit is reasonably successful, and is better than the H- 
dependence derived from Giordano's formulaJ^. It yields 
a critical field H\\ k, 1.1 T. 

To compare with an independent estimate of H\\ based 
on the experimental parameters, we use the dirty limit 
expression for the bulk critical field^^ = (t>o/'K2'^''^\U 
(with ^ the SC coherence length in the dirty limit). Em- 
ploying Eq. (I2H) we observe that, up to the numerical 
constant A/", H\\ is essentially determined by d and ^ only. 
The experimental parameters mentioned in Ref. 'l9'| im- 
ply d ~ 60A, ^ « 1280A (the latter, however, relies on the 
text -book expression for ^ in a dirty Aluminum and the 
measured mean free path). The wires are actually thin 
strips with a rectangular cross section, hence M cannot 
be determined accurately, but is expected to be interme- 
diate between a thin film and a wire, yielding H\\ between 
1.4 T and 2 T. The order of magnitude is consistent with 
the value extracted from Fig. ([3]). Considering the fact 
that the above estimate of H\\ does not rely on any fitting 
parameters, this is a reasonably good agreement. 



IV. CONCLUSIONS 

We have studied the low T resistivity of a dirty SC 
wire using a memory formalism approach. This method 
allows a perturbative treatment of corrections to the low 
energy effective theory, describing the dynamics of phase 



fluctuations in the SC order parameter. In the limits of 
either ideally homogeneous wire or strongly inhomoge- 
neous (granular) wire, our results for p{T) recover the 
power-law behavior ~ T" derived in earlier theoretical 
literature using the instanton technique. We show that 
more generally, the expression for p{T) interpolates be- 
tween the two limits, where the relative weights of each 
are smoothly tuned by a "granularity parameter" D de- 
scribing the typical fluctuations in the diameter of the 
wire. In particular, the inhomogeneity of the wire leads 
to an exponential enhancement of the phase-slips in- 
duced resistivity. 

Following these calculations, we infer that the low en- 
ergy theory is likely to be relevant to experimental mea- 
surements of p{T) in the low T, QPS regime. We di- 
rectly test the validity of the power-law ansatz as a fit to 
experimental data obtained in a variety of different sam- 
ples, and found a very good agreement in sufficiently long 
wires. The exponent a extracted from the p{T) data is 
found to be consistent with the values obtained from an 
independent estimate, based on known parameters of the 
corresponding samples (as much as the information was 
available). In addition, its dependence on a parallel mag- 
netic field is consistent with the theoretical expectation. 
We conclude that the low-T theory for dirty SC wires 
provides a plausible interpretation of data presented in 
the experimental literature, which is better justified than 
previously used phenomenological effective models. Our 
conclusions will hopefully encourage a further, more sys- 
tematic investigation of the comparison between theory 
and experiment. 
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APPENDIX A: DERIVATION OF THE 
RESISTIVITY WITHIN THE MEMORY 
FUNCTION APPROACH 

In this Appendix we review the general aspects of the 
memory matrix approach^^i^iS^iS^j which turns out to 
serve as a useful tool for evaluating transport properties 
in systems where approximate conservation laws can be 
easily identified. As we show below, the calculation of 
d.c. electric resistivity of a SC wire provides a particu- 
larly simple example for the application of this approach, 
which amounts to a straightforward perturbative treat- 
ment of any processes responsible for relaxing the trans- 
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port current. 

At finite but low T, the fixed point Hamiltonian Ho 
[Eq. ([3])] provides the leading contribution to thermo- 
dynamic properties of the SC wire. However, by itself 
it does not give access to transport properties: being a 
translationally invariant intcgrablc model, it possesses an 
infinite number of conserved currents, including in par- 
ticular the electric charge transport current. Since the 
current cannot degrade, the d.c. conductivity is infinite 
even for T > 0. To get the leading non-trivial contribu- 
tion to transport, it is therefore necessary to add the irrel- 
evant corrections, leading to a slow but finite relaxation 
rate of the currents. In our case, the prominent correc- 
tions are the phase-slips terms H^" [Eq. note that 
other irrelevant terms, e.g., of the form J(9(/))"'((96')", 
J {d™'(p){d(f>)^ with n + m > 2 as well as pair break- 
ing terms, which are already neglected in Eq. 12]), do 
not contribute significantly to transpor t^^'^^ . We then 
evaluate the resistivity p{T) employing the memory ma- 
trix approach, which takes advantage of the fact that the 
memory matrix - a matrix of decay-rates of the slowest 
modes in the system - is perturbative in the irrelevant op- 
erators. This is in contrast with the conductivity matrix, 
which is a highly singular function of these perturbations. 

A crucial step in the derivation of transport coefficients 
in the memory matrix approach is the identification of 
primary slowly decaying currents of the system. These 
are conserved "charges" of the fixed point Hamiltonian 
Hq, whose conservation is slightly violated by certain 
irrelevant perturbations. In the case of the Luttinger 
model^^, these include the charge current Je [Eq. ([8])] 
and, in addition, the total translation operator 



Jt = — dxlldx 



(Al) 



In the above, we approximate the currents by contribu- 
tions from the collective degrees of freedom (phase fiuc- 
tuations) only. Non-superconducting contributions as- 
sociated with unpaired electrons are exponentially sup- 
pressed as e-'^^'^ (with A - for T < T^. The correla- 
tors of Je, Jt determine the conductivity matrix a{Lu,T) 
at frequency uj and temperature T via the Kubo formula; 



(c.,T) = 



TL 



dte^'^' {Jp{t)\jg) 



(A2) 



where following Ref. [23| we have introduced the scalar 
product (of any two operators A and B) 



{A{t)\B) 



T 



l/T 



d\{A{tyB{iX)) . (A3) 



The d.c. electrical resistivity is then given by p{T) ~ 
l/aee{T), where fT(T) is the w — > limit of aed^^jT). 
However, a direct application of Eq. (|A2[) at low T, ui is 
rather subtle: since [Jq,Ho] = (for q = e,T), the relax- 
ation rate of the currents dtJg = i[Jg,H] is dictated by 
the irrelevant corrections, hence tends to vanish in the 



limit T — > 0. This leads to divergences in the conductiv- 
ities, since the currents do not decay in this limit. 

To enable a controlled perturbative expansion in the 
relaxation rates dtJq, we therefore recast the conductivity 
matrix in terms of a memory matrix M: 



in which 



Mpg{UJ,T) 



1 

TL 



dtJp 



Q- 



QCQ' 



dtJg] (A5) 



and X is the matrix of static susceptibilities 



1 



(A6) 



Here £ is the Liouville operator defined by CA = [H, A] , 
and Q is the projection operator on the space perpendic- 
ular to the slowly varying variables Jq, 



Q=l-J2\J^)±ix-\qiJ, 



(A7) 



Note that similarly to Ref. [301, we choose a convenient 
definition of the memory matrix [Eq. (|A5[) ]. which is 
slightly different than the standard literature. The per- 
turbative nature of M is transparently reflected by this 
expression: in particular, the operators dtJq are already 
linear in the irrelevant corrections to Hq. This enables 
a systematic perturbative expansion of the correlator in 
Eq. (|A5|) in the small parameter characterizing the rela- 
tive size of the corrections - in our case, the exponential 
factors [see Eq. (g])]. 

We now obtain an expression for the d.c. electric re- 
sistivity by setting = in Eq. (|A4p . This yields 



PiT) 



XieMTT + 2XeTXeeMeT + XtrMe. 



(AS) 



where M is evaluated from Eq. (|A5p in the limit a; — > 0. 
(Note that here we have used the fact that both matrices 
M{T) and x{T) are symmetric.) The relaxation rate 
operators dtJq are given by 



(A9) 



where we have defined the 'force' operators 

^;\„EEz[J„77r] • (AlO) 



To leading order in the perturbations H^'', Eq. (|A5[) 
for M is greatly simplified: one can set C = Co with 
£0 = [Ho, .] and Q = 1. This yields 



(All) 
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where the elements of Mps,n are given by 



hm 



ipp -pi \0 

\ ps,n' ps,n/u;— 



(A12) 



in which is the retarded correlation function 

evaluated with respect to i?o- 

Substituting Eqs. and (gl) in Eq. (lAlOp and tak- 

ing the L ^ oo limit, it is easy to see that „ = 0. 
Indeed, this follows from the translational invariance of 
HP^. As a result, Mtt and MeT identically vanish. This 
implies that i/both XeT and M^e do not vanish, Eq. (|A8[) 
yields /9(T) = even at finite T; namely, the creation of 
free phase-slips appears to be insufficient to generate a 
finite dissipation. Such a result cannot be reconciled with 
our understanding, that dissipation should occur at the 
normal cores during a phase-slips event, where the wire 
behaves temporarily as a normal dirty metal. Indeed, 
we should recall that the translational invariance of H^^ 
is not of microscopic origin; rather, these terms in the 
effective Hamiltonian are obtained after averaging over 
disorder in the electron system^. Their form reflects the 
total absorption of the finite momentum P = 2'Ksns gen- 
erated in a phase-slip process by the core electrons. In 
a clean SC wire, such momentum transfer is not effec- 
tive, leading to an exponential suppression of Mtt, M^e 
and yielding a finite but exponentially small resistivity 
p(T) ~ (see Ref. [TSIllI)- Stronger disorder 

in the underlying electronic system should only enhance 
the resistivity, and not vice versa! This apparent paradox 
is resolved by the remarkable observational that, follow- 
ing the commensurate relation between the total momen- 
tum P and density per unit length, XeT actually vanishes 
identically, even at finite T. Setting XeT = in Eq. (jASp . 
we then find that the expression for the resistivity p{T) 
reduces to the simplified form Eq. ([9]). 

The above arguments are accurate as long as indeed 
the currents Je, Jt are given by Eqs. ([5]), (|Aip . i.e.. 



when normal electron contributions to the currents are 
neglected. However, if such normal contributions are not 
negligible, they would at the same time modify the trans- 
lation invariance of the Hamiltonian, and all entries in 
Eq. (|A8|) would be finite. Similarly, translational in- 
variance is explicitly broken once we account for random 
spatial fluctuations in the SC wire cross section, by in- 
troducing the inhomogeneous phase-slip term Eq. (|2ip . 
This modifies the force operator Fp^ [now given by Eq. 
p4| ]. and induces a non-trivial contribution to F^^: 



F' 



.2 2 

-I— V 



dx y{x) {dj) sin (26*) 



However, the resulting matrix M{T) obtained after sub- 
stitution in the approximate form Eq. (|A12p is still di- 
agonal: to see this, we note that the off-diagonal element 
is given (after disorder averaging) by 



K 



(A13) 



X I dtt I dx e-^(^)e-^(°) lin{CeT{x,t)} 



where 



(A14) 



in which Ge{x,t) is given by Eq. (llSp . It is apparent 
from Eqs. (|A14p and that CeTix,t) is an antisym- 
metric function of x, and hence the integral over x in 
Eq. (jAl3P vanishes, yielding MeT — 0. The static sus- 
ceptibility matrix [Eq. (jA6P ] evaluated to leading order 
in the perturbation H^p is also approximately diagonal, 
with Xee given by Eq. ([M]) . We therefore find that the 
leading contribution to Eq. (jASP for the resistivity again 
practically reduces to Eq. 
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